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Normal  Shock  Interaction  With  a Turbulent 

★ 

Boundary  Layer  on  a Curved  Wall 
6.  R.  Ingerh  and  H.  Sobieczky++ 


Abstract 


A detailed  analysis  is  made  of  longitudinal  surface  curvature  effects  on 
the  Interaction  of  a weak  normal  shock  with  a non-separating  two-dimensional 
turbulent  boundary  layer.  It  is  shown  that  the  Interactive  viscous  displace- 
ment effect  on  the  local  outer  inviscid  transonic  flow  completely  eliminates 
the  well-known  singularity  pertaining  to  purely  Inviscid  flow  on  a curved  wall, 
i.e.,  the  interactive  pressure  field  Is  regular  behind  the  shock.  A study  of 
the  inner  interaction  solution  within  the  boundary  layer,  however,  reveals 
that  curvature  can  influence  the  interaction  but  for  a hitherto-overlooked 
reason:  the  effect  on  the  turbulent  eddy  viscosity,  which  alters  the  boundary 
layer  profile  shape  upon  which  the  interaction  depends.  An  approximate  non- 
asymptotlc  solution  is  given  which  Incorporates  this  effect  and  example  num- 
erical results  are  presented  and  verified  by  comparison  with  experimental 
data.  Small  amounts  of  curvature  (K.fi  -.01-. 02)  are  found  to  moderately  spread 
out  and  thicken  the  Interaction  zone  while  also  beneficially  delaying  the  onset 
of  any  incipient  separation  that  occurs  under  the  shock  foot. 
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Nomenclature 


Sonic  speed  of  sound 

Skin  frlctloi.  coeff.  <■  VA.  U4  > 

Functions  defined  In  Eqs.  4,  7 

Body  shape  function,  see  Fig.  4 

Incompressible  form  factor  (=  ) 

1 1 ' 

Wall  curvature 

Distance  to  shock  location  (Fig.  7) 

Mach  number 
Static  pressure 

Wall  radius  of  curvature  (=  K"1) 

Reynolds  number  [-/>  L/y^  ) 

Transformed  dependent  variables,  Eq.  1 

Velocity  components  In  x,y  directions,  respectively 

Resultant  flow  velocity 

Coordinates  along  and  normal  to  body  surface 

Stretched  independent  variables,  Eqs.  2 and  5 

Shock  wave  angle.  Fig.  3 

Specific  heat  ratio 

Boundary  layer  thickness 

Displacement  thickness 

Non-dimensional  interactive  displacement  thickness.  Fig.  7 
Resultant  flow  direction  angle 
Coefficient  of  viscosity 
Similarity  parameter,  Eq.  2 


p Density 

t Shear  stress 

9*  Momentum  thickness 

Sub  and  Super-scripts 

0 Undisturbed  incoming  boundary  layer  properties 

1 Inviscid  flow  conditions  ahead  of  shock 

2 Inviscid  flow  conditions  behind  shock 

e Conditions  at  boundary  layer  edge 

w Conditions  at  wall  surface 

( )'  denotes  perturbation  from  "0"  state 
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1 . INTRODUCTION 

The  Influence  of  surface  curvature  on  transonic  shock-boundary  layer  Inter- 
action is  an  important  basic  and  practical  question,  since  these  interactions 
often  occur  on  curved  surfaces  and  have  significant  effects  on  the  flow  around 
aerodynamic  bodies  and  cascade/turbine  blades  and  on  internal  flows  within  chan- 
nels and  diffusors.  As  Is  well-known,  Oswatitsch  and  Zierep1  studied  the  related 
problem  of  normal  shock  impingement  on  a curved  wall  in  a purely  inviscid  flow; 
they  found  that  convex  curvature  introduces  a logarithmic  singularity  In  the  wall 

pressure  In  the  form  of  a sharp  post-shock  expansion  (Fig.  1)  and  this  phenomenon 

2 3 

has  indeed  been  observed  in  inviscid  numerical  solutions  by  Emmons  , Jameson 
and  others.  However,  in  real  flows  even  at  high  Reynolds  numbers,  viscous  smear- 
ing and  upstream  influence  effects  arise  due  to  the  thin  boundary  layer  along  the 

4 5 

surface  which  profoundly  influence  the  physics  of  the  local  Interactive  flow  ’ . 
Consequently,  In  view  of  the  well-known  role  often  played  by  viscosity  in  el- 
iminating singularities,  one  would  expect  the  boundary  layer  to  significantly 
alter  and  perhaps  even  eliminate  this  inviscid  curvature-induced  singularity. 

Yet  such  a singularity  continues  to  be  cited  as  an  explanation  of  various 
features  observed  in  transonic  shock-turbulent  boundary  layer  experiments  along 
curved  walls  and  also  has  been  used  (erroneously  In  our  view)  as  an  outer  bound- 
ary condition  in  a proposed  viscous  interaction  theory  for  curved  walls*’;  thus 
a proper  treatment  of  the  role  of  wall  curvature  In  real  viscous  flows  still 
remains  to  be  given. 

In  the  present  paper,  we  address  this  issue  for  the  case  of  a weak  nearly- 
normal  transonic  shock  interacting  with  a 2-D  non-separating  turbulent  boundary 
layer.  Following  a brief  overview  of  the  general  physical  features  of  the 
problem,  we  first  make  a detailed  analysis  of  the  local  mixed  transonic  inviscid 
flow  structure  outside  the  boundary  layer  when  it  is  interactively  coupled  with 
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the  attendant  thickness  response  of  the  layer.  The  results  show  that  the  viscous 
displacement  effect  eliminates  the  aforementioned  wall -curvature  singularity, 
i.e..  It  in  fact  a mathematical  creature  of  a purely  inviscid  model  and  the 
true  influence  of  curvature  of  the  outer  flow  is  actually  quite  well-behaved. 

We  then  examine  the  inner  disturbance  flow  within  the  boundary  layer,  where  it 
is  shown  that  surface  curvature  does  exert  an  influence  on  the  interaction,  but 
from  a hitherto-overlooked  source:  its  effect  on  the  turbulent  eddy  viscosity. 
Using  an  approximate  non-asymptotic  method  of  solution  to  account  for  this  effect, 
example  numerical  results  are  given  for  typical  interaction  and  shown  to  be  in 
accord  with  Ackeret,  Feldman  and  Rott's  experimental  observations^.  A re- inter- 
pretation of  curved  wall  data  in  the  light  of  the  present  findings  is  also  given. 

2.  TYPICAL  FEATURES  OF  THE  INTERACTION 

As  a background  for  the  present  approach,  we  first  briefly  outline  the 
principal  physical  features  of  a non-separating  normal  shock-turbulent  boundary 
layer  interaction.  To  fix  ideas,  we  consider  an  undisturbed  equilibrium  tur- 
bulent  boundary  layer  in  the  Reynolds  number  range  10  < Re^  < 10  that  is 
perturbed  by  a normal  or  nearly-normal  shock  weak  enough  to  avoid  separation 
(M.|  $ 1.3);  whether  the  wall  is  flat  or  curved,  it  is  an  experimental  fact  under 
these  conditions  that  the  resulting  interaction  flow  pattern  is  of  the  relatively- 
simple  type*  illustrated  in  Fig.  2A.  The  boundary  layer  spreads  out  the  shock 
perturbation  upstream  as  well  as  downstream  to  give  a wall  pressure  distribution 
such  as  shown  in  Fig.  2B.  There  is  observed  a corresponding  significant  local 


* When  separation  occurs,  this  pattern  changes  rapidlv  to  a much  more  com- 
plicated one  involving  a bifurcated  shock  structure^. 


growth  of  the  boundary  layer  displacement  thickness,  particularly  behind  the  shock 
(Fig.  2C) , so  that  even  though  the  original  undisturbed  boundary  layer  is  thin  and 
negligibly  affects  the  overlying  inviscid  flow,  this  is  no  longer  the  case  in  the 
interaction  zone  astride  the  shock  whether  or  not  the  wall  is  flat. 


Thus,  the  local  outer  inviscid  flow  "sees"  the  overall  viscous  effects  in 
the  interaction  as  an  effective  thickening  of  the  geometrical  surface  by  an  amount 
equal  to  the  Interactive  displacement  thickness  growth.  In  the  vicinity  of  the 
shock,  this  is  seen  to  involve  (a)  a rapid  growth  in  the  effective  body  slope  at 
the  shock  foot  followed  by  (b)  an  overall  downstream  Increase  of  the  effective 
body  thickness.  At  the  very  least,  these  two  features  must  be  included  in  any 
realistic  mathematical  model  of  the  true  outer  inviscid  flow  in  such  interactions. 

3.  OUTER  INVISCID  FLOW  MODEL 

Consider  the  local  mixed  inviscid  transonic  flow  in  the  neighborhood  of  the 
shock  outside  the  boundary  layer,  assuming  the  shock  weak  enough  to  be  isentropic. 
Then  If  W and  are  the  magnitude  and  direction  angle,  respectively,  of  the 
resultant  flow  velocity  and  we  introduce  the  transformed  dependent  variables 


o ■ / 2 i / u 

S = + 2 (y+1 ) • o ^ -*  - U 

t = o'1  U 


as  functions  of  the  stretched  independent  variables 

X = x (t 

Y = [f  y ^ 

where  the  constant  T is  chosen  later  as  an  appropriate  similarity  parameter, 

then  the  governing  transonic  flow  equations  are: 

cl/3  3S  3t  _ n (■ 

53  3X  ” 3Y  " u 
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This  pair  of  equations  in  general  has  either  supersonic  wave-type  solutions 
(S  > 0,  minus  sign  in  3b)  or  subsonic  Cauchy  - Riemann  - type  solutions  (S  < 0, 
positive  sign).  If  we  further  restrict  attention  to  small  perturbations  about 
the  shock  conditions,  Eqs.  (3)  become  linear  in  the  leading  approximation. 

We  now  choose  (T  such  that  upstream  of  the  shock  S = S1  =1;  then  linear- 
ization of  Eq.  (3)  gives  a wave  equation  with  solutions  of  the  forn/ 


S = S1  + F] (c)  + Gj (n) 


t = t1  + F^c) 


G-j  (n) 


(4) 


where  c = X + Y and  n = X - Y.  We  linearize  the  subsonic  equations  for  per- 
turbations of  the  post-shock  condition  S2  < 0,  t 2 and  find  the  solution  of  the 
Cauchy-Riemann  equations^: 


Z = X + iS21/3Y 
= S2  + it2  + F2(Z) 


(5) 


A particular  case  of  a pair  of  local  solutions  (4)  and  (5)  pertains  to  the  flow 
ahead  and  behind  an  oblique  shock  of  the  "strong  branch"  type,  i.e.,  with  sub- 
sonic flow  downstream.  The  transonic  shock  jump  conditions  link  the  pre-  and 
post-shock  conditions  as  follows: 


St  = t2  - t,  = l'm  . 3 . (1  - S22/3)'/2  0 * S22/3)  (6a) 

Ctg  6 • 2-5'6  . 31'3  (y  ♦ 1)'/3  «’/3  (1  - S22'3>'/2  (6b) 

★ 

Thus,  if  at  some  point  P there  is  a supersonic  flow  W-j  > a at  inclination 
V-j  that  is  decelerated  by  an  oblique  shock  to  some  W2  < a at  angle  (Fig. 
3A)  then  W1 , *1^,  W2,\£  are  related  in  terms  of  Sy  ty  Sy  1 2 by  means  of 


Eqs.  (la,  lb,  6a).  This  is  the  transonic  shock  polar  (Fig.  3b)  and  we  find 

★ 

the  wave  angle  B corresponding  to  the  flow  deflection  angle  from  Eq.  (6b). 

Then  with  given  shock  jump  conditions  at  P,  the  local  solutions  (4)  and  (5) 

may  be  used  u determine  the  velocity  gradients  as  well  as  the  curvature  of 

the  shock  shape  in  the  vicinity  of  P. 

This  local  transonic  inviscid  perturbation  field  now  can  be  linked  to  the 

inner  viscous  disturbance  flow  by  noting  from  above  that  the  outer  flow  sees 

2 

an  effective  body  consisting  of  the  bare  wall  hw  = h^  = C^X  plus  a slope  jump 
C2  at  the  shock  foot  followed  by  a further  interactive  displacement  thickness 
growth,  giving  h2  = hw  = C^X  + D^X  (see  Fig.  4a).  Upstream  of  the  shock, 
an  arbitrary  linear  variation  of  inviscid  flow  Mach  number  along  the  wall  is 
assumed  (Fig.  4b)  to  allow  for  either  decelerating,  accelerating  or  uniform  in- 
coming supersonic  flow.  Now  in  terms  of  this  flow  model,  a local  analytic 

7 ★★ 

solution  can  be  carried  out  by  a double  series  expansion  analysis  of  the  up- 
stream and  downstream  flow  properties  about  the  shock-foot  origin  using  the 
following  particular  solutions  of  Eqs.  (4)  and  (5): 

F,  ■ A,s 

G1  * ** Bln  (7) 

F2  = A2Z  + B2Z  In  Z 

The  downstream  (complex)  coefficients  A2>  B2  are  related  to  their  upstream 
(real)  counterparts  and  the  shock  shape  parameters  by  the  transonic  shock  polar. 
The  result  of  this  analysis  yields  two  basically  different  types  of  physical 

* It  is  assumed  here  that  the  deflection  angle  -}£  lies  below  the  maxi- 
mum value  of  the  shock  polar,  which  is  reasonable^in  practice  if  the  in- 
cident shock  is  weak  enough  to  avoid  separation. 

**  The  realistic  assumption  is  made  here  that  the  radius  of  any  wall  curvature 
is  large  compared  with  a typical  boundary  layer  displacement  thickness. 
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behavior  (one  regular,  the  other  singular)  depending  on  whether  the  interactive 
deflection  effect  is  included  (C2  >0,  f ) or  neglected  a priori  (C^  = 0, 

D2  = Cj ) ,as  follows. 

(1)  C?  > 0,  d2  f c1 

Regular  Solution  F2:  f 0,  B2  = 0 (8a) 

Parabolic  Oblique  Shock  Shape:  = Qy  + Py^  (8b) 

Linear  Downstream  Wall  Pressure:  M2~l  = K2  + L^X  (8c) 

as  typically  illustrated  in  Fig.  4b  (expressions  for  P,  Q,  K2>  L2  are  given 

in  trie  Appendix).  Owing  to  the  displacement  thickness  effect,  note  here  that 

the  inviscid  normal  shock  actually  is  curved  near  the  boundary  layer  edge  with 

g 

a slightly  oblique  angle  at  the  foot. 

(2)  C2  = 0,  D2  = C1 


Singular  Solution  F2>- 

a2  = o,  b2  f 0 

(9a) 

Curved  Shock  Normal  to  Wall: 

xsh  = y2  (Q  + R ln  y) 

(9b) 

Logarithmically  Singular  Pressure: 

M2=-K]  (1-E2X-F2X  ln  X) 

(9c) 

as  schematically  depicted  in  Fig.  5 (the  appropriate  coefficients  can  be 
found  in  Ref.  8.)  This  "bare"  curved  wall  solution  is  in  fact  the  Oswatitsch- 
Zierep  singular  solution^  mentioned  earlier.  It  has  been  of  use  mainly  to  verify 
computational  results  for  inviscid  flows  with  recompression  shocks:  fully  con- 
servative codes  usually  give  results  with  a post-shock  expansion,  demcs. rating 
the  locally-exact  treatment  of  the  recompression  shock  if  the  computational  grid 

3 

is  fine  enough  . 

The  authors  of  Ref.  1 were  inspired  to  clarify  the  problem  of  the  normal 

4 

shock  on  a curved  wall  by  Ackeret's  remarks  about  an  experimentally-observed 


6 


sharp  post  shock  expansions  at  the  boundary  layer  edge.  But  these  were  observed 

for  luminar  boundary  layer  interactions  where  a lambda  shock  configuration  occurs 

and  an  expansion  fan  follows  the  strong  main  shock  sitting  on  a large  separation 

bubble  (Fig.  6).  For  the  experiments  with  turbulent  boundary  layer  inter- 

4 

action  without  separation,  however,  this  does  not  occur;  instead,  an  oblique 
shock  shape  occurs  on  a ramp-like  thickening  of  the  turbulent  boundary  layer7. 
From  these  observations  we  conclude  that  the  singular  solution  involving  the  log- 
arithmic terms  (Eqs.  9)  might  be  applicable  in  modeling  the  outer  flow  near  the 
shock  foot  in  the  case  of  a fully-separated  flow  where  the  "wall"  shape  is  given 
by  the  curved  separation  bubble  geometry  and  not  by  the  underlying  solid  wall 
(although  in  this  situation,  shock  reflection  is  governed  by  a free  streamline 
instead  of  a solid  wall  boundary  condition  so  that  a singularity  still  need  not 
necessarily  occur).  However,  for  non-separating  turbulent  boundary  layer  inter- 
actions, the  regular  solution  (8)  is  the  appropriate  inviscid  flow  model  even 
on  a curved  wall  because  it  takes  into  account  the  flow  deflection  and  shock 
distortion  due  to  the  interaction-induced  boundary  layer  thickening.  Thus  the 
turbulent  interaction  outer  flow  model  is  primarily  determined  by  the  deflection 
ramp  and  not  by  wall  curvature;  curvature  does  have  some  influence  (as  we  show 
below)  but  it  is  regular  and  does  not  derive  from  any  singular  be- 

havior in  the  inviscid  part  of  the  flow.  This  conclusion  is  concordant  with 
all  known  theory  and  experiment  on  viscous-inviscid  interactions  in  general: 
coupling  of  the  inviscid  flow  with  the  interactive  displacement  effect  in- 
variably eliminates  local  mathematical  singularities  that  otherwise  may  occur 
in  solving  the  flow  equations. 

The  foregoing  discussion  shows  that  Eqs.  (8)  are  the  correct  mathematical 
model  of  the  local  interacting  inviscid  flow  near  the  shock  including  any 
surface  curvature.  As  indicated  in  the  Appendix,  the  downstream  solution 
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coefficients  C2>  D2,  K^,  L2  are  related  to  the  upstream  coefficients  C-j , K1 , 

L1  and  those  of  the  shock  shape  (P,Q)  by  four  equations;  since  the  basic  wall 
shape  (C^)  and  incoming  supersonic  flow  (K^ , L^)  are  usually  known  in  practice, 
any  pair  of  the  remaining  coefficients  must  be  given  to  determine  the  others. 
This  implies  three  possibilities:  (a)  a prescribed  shock  shape  (P,Q)  based 
on  a Schlieren  photo,  leading  to  the  corresponding  effective  interactive  wall 
shape  and  downstream  flow;  (b)  a given  post-shock  Mach  number  distribution 
(«2,  L2)  obtained  from  wall  pressure  data,  leading  to  the  corresponding  shock 
and  effective  wall  shapes;  (c)  a known  post-shock  displacement  effect  (C2>  D2) 
given  either  by  measurement  or  an  interactive  boundary  layer  theory,  leading 
to  the  outer  shock  shape  and  invisc’d  pressure  distribution.  Here  we  are  con- 
cerned  with  (c);  further  discussion  of  the  others  is  given  elsewhere  ’ . 

4.  INNER  VISCOUS  FLOW  MODEL 


4.1  General  Features 

The  present  treatment  of  the  interactive  perturbation  field  within  the 
inner  turbulent  boundary  layer  region  is  based  on  an  appropriate  extension  of 
Inger  and  Mason's  approximate  analytical  theory^  to  include  consideration  of 
wall  curvature  effects.  Consider  a known  adiabatic  boundary  layer  profile  MQ(y) 
subjected  to  small  transonic  disturbances  due  to  an  impinging  weak  and  nearly 
normal  shock.  In  the  practical  Reynolds  number  range  of  interest  here 

c o 

(10  £ Re^  <_  10  ),  we  employ  a non-asymptotic  disturbance  flow  model  in  the 

turbulent  boundary  layer  patterned  after  the  Lighthill-Stratford-Townsend  double- 
23  24 

deck  approach  ’ that  has  proven  highly  successful  in  treating  a variety  of 
other  problems  involving  turbulent  boundary  layer  response  to  strong  rapid  adverse 
pressure  gradients.  This  was  done  because  of  the  large  body  of  turbulent  bound- 
ary layer-shock  interaction  data  that  strongly  supports  such  a model  in  this 
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25  26 

Reynolds  number  range  * (including  the  present  transonic  regime)  and  the 

27 

findings  of  a separate  general  theoretical  study  suggesting  that  the  results 

20  28 

of  the  asymptotic  theory  approach  for  very  high  Reynolds  numbers,  ’ although 
rigorous  in  this  limit,  may  not  extrapolate  down  to  the  present  Reynolds  number 
range. 

The  resulting  flow  model  consists  of  an  inviscid  boundary  value  problem 
surrounding  a nonlinear  shock  discontinuity  and  underlaid  by  a thin  viscous  dis- 
turbance sublayer  as  schematically  illustrated  in  Fig.  7.  An  approximate  ana- 
lytical solution  is  further  achieved  in  the  leading  approximation  by  introducing 
some  engineering  simplifications^  (including  the  assumption  of  small  linearized 

disturbances  ahead  of  and  behind  the  nonlinear  shock  jump  plus  neglect  of  the  de- 

12 

tailed  shock  structure  within  the  boundary  layer  ) which  have  been  shown  to  give 

accurate  results  for  the  overall  properties  of  engineering  interest  provided 

★ 

is  not  too  close  to  unity  (M^  > 1.05  ).  This  yields  a triple-deck  solution  con- 
sisting of  (a)  linearized  potential  supersonic  disturbance  flow  in  region  1 plus 
the  subsonic  disturbance  flow  in  region  3 caused  by  the  interaction-generated 
interface  displacement  ^(x)  anc*  the  post-shock  perturbations  along  x = 0+  due 
to  the  impongement  of  region  1 Mach  wave  disturbances  on  the  shock;  (b)  under- 
neath, a rotational  inviscid  disturbance  flow  region  2 in  the  boundary  layer;  (c) 

1 3 

a thin  viscous  disturbance  sublayer  within  the  linear  part  of  the  basic  velocity 

profile  U = (t  /m  )y  which  carries  most  of  the  upstream  influence  of  the  inter- 

1 3 

action.  A known  solution  for  this  sublayer  can  be  extended  to  the  present 


5 8 

* As  far  as  the  overall  interaction  soi.'cion  for  10  < ReL  < 10  is  concerned, 

the  nonlinear  shock  jump  conditions  plus  the  non-uniform  viscous  flow  effects 
in  the  boundary  layer  reduce  the  lower  Mach  number  limit  otherwise  pertaining 
to  the  linearized  supersonic  theory  in  purely  inviscid  potential  flow  (see  Ref. 
21  for  more  detailed  discussion). 


« 
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problem  to  obtain  the  displacement  effect  (effective  wall  position  seen  by  the 
overlying  rotational  inviscid  disturbance  flow)  and  also  the  corresponding  dis- 
turbance skin  friction  caused  by  the  pressure  disturbance  field. 

The  matching  of  these  regional  solutions  yields  linear  integral  equations 
that  can  be  solved  by  operational  methods  for  the  disturbance  pressure  along 
both  the  boundary  layer  edge  and  wall^.  The  remaining  Interactive  flow  prop- 
erties can  then  be  determined  in  terms  of  p^,  including  the  interactive  growth 
in  displacement  thickness  and  the  corresponding  skin  friction  perturbation.  It 
is  noted  that  this  latter  perturbation  solution  has  been  extended  downstream  as 
well  as  upstream  of  the  shock  and  further  corrected  for  the  nonlinear  inertia 
effects  in  an  adverse  pressure  gradient  using  the  general  non-dimensional  wall 

5 

shear-pressure  solution  ahead  of  separation  given  by  formal  triple  deck  theory  ; 

when  converted  to  turbulent  flow  by  expressing  all  results  in  terms  of  in- 

o 

stead  of  ReL  , this  yields 

C'  F(x/6  ) (10) 

*w 

where  the  non-dimensional  function  F is  essentially  unity  ahead  of  the  shock 

-1/3 

x < 0 and  vanishes  slowly  behind  it  with  F ~ (x/6Q)  far  downstream.  It  is 

seen  from  Eq.  (10)  that  depending  on  Reynolds  number  (Cf  ),  a sufficiently  strong 

o 

interactive  pressure  rise  can  cause  incipient  separation  (C^  -*■  0)  near  the  shock. 

This  flow  model  contains  all  the  essential  global  features  of  the  mixed 

transonic  character  of  the  non-separating  normal  shock-turbulent  boundary  layer 

interaction  problem  including  the  significant  lateral  pressure  gradient  effects, 

upstream  influence  and  interactive  skin  friction  for  an  arbitrary  input  turbulent 

14  15 

boundary  layer  profile.  Detailed  comparisons  with  experiment  * have  shown 
that  it  gives  a very  good  account  of  all  the  important  engineering  features  of 


10 


I 


the  Interaction  over  a wide  range  of  Mach-Reynolds  number  conditions  at  very  low 
computational  cost.  Hence  the  theory  provides  a sound  basis  for  interpreting 
experimental  data  on  unseparated  flow  and  for  further  extensions  to  account  for 
slight  wall  curvature. 

It  is  noted  that  this  solution  in  principle  may  be  used  with  a mean  turbulent 
boundary  layer  profile  input  from  either  experimental  measurement  or  any  theoret- 
ical prediction  method.  Here  we  have  chosen  an  approximate  but  accurate  and  par- 
ticularly convenient  analytical  profile  model  by  LeFoll  that  is  in  frequent  use 
in  boundary  layer  studies,  has  considerable  parametric  generality  and  proves 
especially  well-suited  to  perturbation  problems  in  non-uniform  flows  (a  detailed 
description  may  be  found  elsewhere16,  17).  It  is  characterized  quite  generally 
by  three  parameters  (M^ , boundary  layer  thickness  Reynolds  number  and  the  in- 
coming shape  factor)  and  readily  permits  the  introduction  of  the  effects  of  wall 
curvature  as  shown  below. 

4.2  Longitudinal  Curvature  Effects 

Now  wall  curvature  influences  the  inner  solution  in  two  ways:  (1)  it  alters 
the  incoming  undisturbed  turbulent  boundary  layer  profile  upon  which  the  per- 
turbation solution  depends;  (2)  it  introduces  new  explicit  terms  in  the  small 
disturbance  equations.  The  first  is  by  far  the  most  important  and  is  discussed 
further  below;  the  secondary  effects  (2)  are  in  fact  negligible  for  the  small 
relative  curvatures  U/RDj£.02)  encountered  in  practice  according  to  the 
following  considerations.  As  regards  the  (regular)  influence  of  curvature  on  the 
inviscid  outer  boundary  conditions  via  the  displacement  effect,  the  foregoing 
analysis  shows  that  this  is  of  order  6Q*/RB  and  hence  very  small  in  the  leading 
approximation.  Likewise,  under  the  continued  assumptions  that  the  interactive 
disturbances  remain  uncorrelated  with  the  background  turbulent  fluctuation  field 
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and  that  the  viscous  disturbance  sublayer  lies  within  the  laminar  sublayer 

region  of  the  turbulent  boundary  layer,  the  new  terms  in  both  the  rotational 

inviscid  and  viscous  disturbance  sublayer  equations  can  be  shown1 ^ to  be  of 

the  same  order  or  smaller  and  hence  also  negligably  small.  Thus,  to  a con- 

sistant  degree  of  first  approximation,  the  form  of  the  equations  governing 

the  disturbance  flow  in  the  boundary  layer  is  not  altered  by  wall  curvature. 

Regarding  the  basic  turbulent  flow  pro-file,  once  again  the  explicit  new 

curvature  terms  in  the  flow  equations  governing  this  profile  (including  the 

centrifugal  lateral  pressure  gradient  effect)  have  all  been  shown  to  be  of  the 

order  6Q/Rg  and  hence  negligible  (e.g.,  see  Bradshaw's  monograph  and  more  re- 
19 

cent  work  ).  However,  the  influence  on  the  eddy  viscosity  relation  is  known 
to  have  an  order  of  magnitude  larger  effect  (10  to  20  * 6 /RD ) on  the  overall 

0 D 

skin  friction t and  form  factor  H’(Ref.  18);  since  it  has  been  firmly  established 
o 

that  the  role  of  the  profile  in  the  interaction  mainly  derives  from  these  proper- 
ties1^, this  influence  is  therefore  the  cause  of  any  significant  curvature  effect 
on  the  interaction  solution. 

By  virtue  of  the  foregoing  arguments,  the  curvature  effect  can  be  reasonably 
estimated  by  incorporating  into  the  incoming  boundary  layer  profile  model  ap- 
propriate corrections  for  the  influence  of  curvature  on  C/0,  the  form  factor  and 
thickness.  Examination  of  Bradshaw's  comprehensive  study  shows  that  over  a 
wide  range  of  parametric  and  local  conditions,  these  corrections  for  small  long- 
itudinal curvature  and  non-separating  flows  can  be  represented  by  the  engineering 
approximations* 


* It  is  emphasized  that  the  use  of  a more  detailed,  sophisticated  model  will 
not  alter  the  general  trends  and  qualitative  conclusions  of  primary  concern 
here. 
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(lib) 


where  to  this  first  order  accuracy  the  corresponding  effect  on  6Q  is  much 
smaller  and  therefore  neglected.  Note  that  the  typical  value  d /R„  % .01 

0 D 


yields  a reduction  and  increase  in  Cf  and  H.  of 

To  71 


10%  and  5%,  respec- 


tively, which  are  an  order  of  magnitude  larger  than  the  explicit  new  curvature 
terms  in  either  the  mean  or  perturbation  flow  equations  (including  the  outer 
inviscid  disturbances).  The  use  of  Eqs.  (11)  with  dQ/RB  as  an  additional  input 
parameter  in  the  Inger-Mason-Panaras  theory^  *^5 * 7  enables  a straightforward  ap- 
praisal of  the  first  order  curvature  effects  on  the  non-separating  normal  shock- 
turbulent  boundary  layer  interaction  problem.  The  results  of  such  a study  will 
now  be  presented. 


5.  DISCUSSION  OF  RESULTS 


5.1  Parametric  Studv 


The  typical  influence  of  curvature  on  the  interactive  wall  pressure  dis- 
tribution is  illustrated  in  Fig.  8;  the  corresponding  displacement  thickness 
growth  and  local  skin  friction  behavior  are  shown  in  Figs.  9 and  10.  Fig.  8 

shows  that  the  curvature  effect  from  has  a very  small  influence  (this  in 

o 

fact  was  found  to  be  true  over  a wide  range  of  - Re^  conditions),  whereas 
that  from  the  form  factor  is  of  controlling  importance.  Since  curvature  in- 
creases H and  hence  reduces  the  incoming  profile  "fullness",  it  acts  to  spread 
out  the  wall  pressure  disturbance,^7  increase  the  upstream, influence,  reduce  the 
downstream  pressure  level  and  (Fig.  9)  thicken  the  downstream  boundary  layer. 

The  corresponding  influence  on  skin  friction  (Fig.  10)  is  especially  significant: 
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although  curvature  reduces  the  upstream  level,  it  increases  the  local  near 
the  shock  foot  owing  to  the  reduced  adverse  pressure  gradient.  Thus  wall  curva- 
ture retards  the  onset  of  incipient  separation  that  tends  to  occur  slightly  down- 
★ 

stream  of  the  shock  . 

The  curvature  effect  on  the  pressure  disturbance  along  the  outer  edge  of 
the  boundary  layer  (Fig.  11)  is  also  of  interest  and  brings  out  several 
fundamental  points.  First,  a small  subsonic  expansion  region  Is  predicted  to 
occur  right  behind  the  shock  regardless  of  the  wall  curvature;  this  was  shown  to 
be  an  inherent  feature  of  the  mixed  transonic  viscous  interaction  flow  along  a 
flat  surface  by  Inger  and  Mason ^ and  is  further  confirmed  by  detailed  numerical 
solutions^.  Thus,  contrary  to  what  is  sometimes  alleged,  wall  curvature  per  se 
is  not  the  cause  of  such  expansion  regions  in  non-separating  turbulent  interact- 
ions. Second,  when  horizontal  distance  is  properly  rescaled  so  that  x~0  (6Q) 
as  is  done  in  Fig.  11,  these  regions  are  perfectly  regular  as  predicted  by  Eqs. 

(8).  Third,  although  not  the  cause,  convex  wall  curvature  does  strengthen  the 
expansion  and  reduce  the  preceding  local  shock  pressure  jump.  We  re-emphasize, 
however,  that  this  is  a result  of  the  curvature  effect  on  the  boundary  layer  eddy 
viscosity  and  has  nothing  whatsoever  to  do  with  any  invlscid  curvature  singularity. 

5.2  Comparison  with  Experiment 

4 

Ackeret,  Feldman  and  Rott's  famed  experimental  study  of  shock-boundary 
layer  interaction  on  a plate  and  wall  in  the  choked  transonic  flow  of  a slightly- 
curved  wind  tunnel  nozzle  provides  some  examples  of  unseparated  turbulent  flow 

* A1 though  the  skin  friction  solution  of  the  present  theory  is  no  longer 
valid  for  separated  flow  [Cf(x)  < 0 over  some  portion  of  the  wall]  it  is 
still  useful  to  Indicate  trends  toward  this  situation,  l.e.,  where  and  when 
incipient  separation  (C,  -*•  0 at  some  x)  first  occurs. 


suitable  for  comparison  with  the  present  theory.  We  have  chosen  those  for  which 


both  wall  and  inviscid  pressure  distributions,  as  well  as  displacement  thickness^ 
are  given. 

It  is  noted  that  direct  comparisons  with  these  data  involve  numerous  un- 
certainties in  converting  to  the  theoretical  variables  (or  vice-versa),  (a) 
the  boundary  layer  thickness  and  hence  the  inviscid  flow  edge  is  only  approxi- 
mately defined;  (b)  the  shock  wave  location  and  shape  are  uncertain  to  within 
.25  to  .50  « 6Q ; (c)  reading  the  curves  where  they  change  rapidly  introduces  in- 
herent error;  (d)  a significant  background  inviscid  pressure  gradient  is  present 
which  beclouds  interpretation  of  the  outer  fringes  of  the  interaction  ?one  and 
the  character  of  the  incoming  "undisturbed"  turbulent  layer  profile;  (e)  the  up- 
stream boundary  layer  history  is  only  partially  understood,  especially  following 
forced  transition  cases,  and  this  together  with  (d)  cannot  be  fully  accounted 
for  in  the  local  interaction  model;  (f)  a significant  channel  blockage  effect 
occurs  from  the  interactive  boundary  layer  thickening,  which  reduces  the  effective 
theoretical  shock  strength  and  hence  the  downstream  interaction  pressure  level 
(this  was  first  identified  independently  by  Inger  and  Panaras1^’  ^7  and  Melnick 
and  Grossman  who  devised  different  but  equivalent  correction  methods  for  it). 

4 

A typical  non-separating  interactive  pressure  field  measured  by  AF  & R 
is  illustrated  in  Fig.  12a;  experimental  pressure  distributions  along  both 
the  surface  and  the  approximate  boundary  layer  edge  are  compared  in  Fig.  12b 
with  the  theoretical  prediction  (corrected  for  the  estimated  interactive  block- 
age effect  using  the  Inger-Panaras  method)^7,  while  the  corresponding  displace- 
ment thicknesses  are  shown  in  Fig.  12c.  We  note  that  the  relevant  non- 
dimensional  curvature  parameter  for  these  experiments  was  rather  small 
(K- • .0063)  so  that  the  predicted  curvature  effect  on  the  interaction  is 

only  slight,  as  indicated.  In  view  of  the  aforementioned  uncertainties  in 
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interpreting  the  data  plus  the  limitations  of  the  theory  (e.g.,  the  11  dis- 
turbance approximations  and  neglect  of  the  detailed  nonlinear  wave  diffraction 
effects  in  the  supersonic  portion  of  the  boundary  laye»&  the  overall  agreement 
is  seen  to  be  quite  good.  In  particular,  the  following  definitive  features  of 
the  interaction  predicted  by  the  theory  are  well -corroborated:  (1)  the  mag- 
nitude, sign  and  streamwise  extent  of  the  lateral  pressure  gradient  effect  both 
ahead  and  behind  the  shock;  (2)  the  existence  of  a long  slow  interactive  pressure 
rise  (algebraic  rather  than  exponential)  downstream  of  the  shock;  (3)  the  over- 
all streamwise  scale  and  upstream  influence  distance;  (4)  the  magnitude  and 
shape  of  the  interactive  displacement  thickness  growth;  (5)  the  local  inviscid 
pressure  jump  across  the  shock  at  the  boundary  layer  edge;  (6)  a non-singular 
inviscid  subsonic  expansion  region  behind  the  shock  due  to  the  viscous-inviscid 
interaction  and  not  surface  curvature  (note  that  the  zero  curvature  theory  in 
this  region  is  actually  closer  to  the  data),  consistent  with  the  regular  inviscid 
post-shock  model . 

From  foregoing  results  it  is  concluded  that  the  inviscid  regions  of  sharp 
pressure  rise  followed  by  expansion  typified  by  Fig.  12a  should  not  be  physically 
interpreted  as  a symptom  of  an  inviscid  curvature  singularity  but  rather  as  in- 
herent features  of  the  viscous  mixed  transonic  interaction  that  occur  even  on 
a flat  surface  and  which  are  only  secondarily-influenced  by  curvature. 

6.  CONCLUDING  REMARKS 

This  investigation  has  sought  to  delineate  the  essential  influence  of 
surface  curvature  on  non-separating  turbulent  boundary  layer- transonic  shock 
interactions,  to  place  in  true  perspective  the  inviscid  curvature  singularity, 
and  to  devise  an  approximate  analysis  including  curvature  that  enables  both 
parametric  study  and  comparisons  with  experiment.  The  resulting  theory  is  in 


fact  tatamount  to  a non-asymptotic  triple  deck  model  of  the  interaction  per- 
taining  to  the  Reynolds  number  range  10  <_  Re^  <_  10  . Although  it  may  there- 
fore be  improved  upon  by  a more  rigorously-detailed  treatment  the  primary  re- 
sults should  remain  the  same,  including  the  following  overall  qualitative 
conclusion:  when  the  viscous  displacement  effect  on  the  inviscid  flow  is 
accounted  for  (as  it  must  by  very  definition  of  an  interaction  problem),  the 
surface  curvature  effect  on  non-separating  interactions  is  non-singular  and 
derives  from  its  influence  on  the  eddy  viscosity  within  the  boundary  layer. 


APPENDIX 


Determination  Of  Solution  Coefficients  For 
Transonic  Flow  Past  A Ramp  Streamline 
With  Compression  Shock  (Ref.  24) 

We  refer  to  Eqs.  (6)  - (8)  plus  Figs.  3 and  4;  it  is  assumed  that  the 
following  coefficients  are  known:  C^  (solid  wall  curvature),  K^(pre-shock 
Mach  number),  L-j  (pre-shock  acceleration),  C2  (post-shock  streamline  inclina- 
tion at  boundary  layer  edge)  and  D2  (post-shock  streamline  curvature  at 
boundary  layer  edge).  Then  in  terms  of  the  auxilary  formulae 


w = 3L1/(2K1 ) 

(A-l) 

v - 3 • r3'2 

(y  + 1)  Kj~3/2  . C1 

(A-2) 

oj  = 3 • 2’3/2 

(y  + 1)  K^372  • C2 

(A-3) 

At  = 2'5/2  • 3 

• (y  + 1)  C2  K]  3/2 

(A-4) 

inversion  of  the  transonic  shock  polar  (Fig.  3B)  on  the  strong  shock  branch 
(“1  £ S2  £ -3'^)  gives 

S2  = - [(1  - 3"3/2)  yj  1 + a • At2  + b*At4  + c-At6  + 3-3/2]  (A-5) 

where  for  y = 1.40,  a = -0.8255,  b = 0.0,  c = 0.0425.  Further,  the  desired 
unknown  coefficients  take  the  following  values: 

P = 2^2  K^2  (1  - S22/3)1/2  [compression  shock  inclination]  (A-6) 

Q = £ • A • «i  [compression  shock  curvature]  (A-7 ) 

K2=  -S/'3  • K-j  [post-shock  Mach  number]  (A-8) 

L2=  2 p K2/  3S2  [post-shock  acceleration]  (A-9) 


:: 


where 


t = uf[4(l  + C2)  + S22/3  + 1]  + (8  C2  + S22/3  + 1) 
p = (c<o  - T)/n 

X = [c ( u - p • S2  ^3)  + (u  + v]/  3(1  - S^3) 


(A-6) 

(A-7) 

(A-8) 

(A-9) 

(A-10) 

(A-ll) 
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Fig.  2 Typical  Features  of  Normal  Shock-Turbulent  Boundary  Layer  Inter 
actions  Without  Separation 
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Curvature  Effect  on  Interaction  Pressure  Along  the  Wall 
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Fig.  10  Curvature  Effect  on  Interactive  Skin  Friction  Distribution 


Fig.  11  Curvature  Effect  on  Interaction  Pressure  Distribution  Along  Boundary 
Layer  Edge 


Fig.  12B  Comparison  of  Present  Theory  with  Experiment-Pressure  Field 
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influence  the  interaction  but  for  a hitherto-overlooked  reason:  the 
effect  on  the  turbulent  eddy  viscosity,  which  alters  the  boundary  layer 
profile  shape  upon  which  the  interaction  depends.  An  approximate  non- 
asymptotic  solution  is  given  which  incorporates  this  effect  and  example 
numerical  results  are  presented  and  verified  hy  rnm|Tari<;nn  with  experi- 
mental data.  Small  amounts  of  curvature^TK  <s  «.01  -To2)^are  found  to 
moderately  spread  out  and  thicken  the  interaction  zone  while  also  benefi- 
cially delaying  the  onset  of  any  incipient  separation  that  occurs  under 
the  shock  foot.j^ 
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